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Photons produced in the annihilations of dark matter particles can be detected by gamma-ray 
telescopes; this technique of indirect detection serves as a cornerstone of the upcoming assault on 
the dark matter paradigm. The main obstacle to the extraction of information about dark matter 
from the annihilation photons is the presence of large and uncertain gamma-ray backgrounds. We 
present a new technique for using gamma-ray data to constrain the properties of dark matter that 
makes minimal assumptions about the dark matter and the backgrounds. The technique relies on 
two properties of the expected signal from annihilations of the smooth dark matter component in 
our galaxy: 1) it is approximately rotationally symmetric around the axis connecting us to the 
Galactic Center, and 2) variations from the mean signal are uncorrelated from one pixel to the next. 
We apply this technique to recent data from the Fermi telescope to generate constraints on the dark 
matter mass and cross section for a variety of annihilation channels. We quantify the uncertainty 
introduced into our constraints by uncertainties in the halo profile and by the possibility that the 
halo is triaxial. The resultant constraint, the flux F < 4.5 x 10 -6 cm -2 s _1 sr" 1 for energies between 
1 and 100 GeV at an angle 15° away from the Galactic Center, translates into an upper limit on the 
velocity weighted annihilation cross section of order 10 -25 cm 3 s _1 depending on the annihilation 
mode. 

PACS numbers: 95.35.+d; 95.85.Pw 

I. INTRODUCTION 

Evidence for the existence of non-baryonic dark matter has been accumulating for many decades. Combined 
constraints from measurements of anisotropics in the cosmic microwave background radiation, the shape of the galaxy 
power spectrum, and the Hubble constant fix both the total matter and the baryon densities, indicating that non- 
baryonic dark matter makes up 85% of the matter density of the universe Despite the preponderance of evidence 
for its existence, however, little is known about the identity of the dark matter. One way to glean information about 
the properties of this mysterious substance is through indirect detection, the observation of the annihilation products 
of dark matter particles. 

Indirect detection is an attractive prospect for several reasons. Its primary advantage is that, while dark matter 
itself is very difficult to detect, the annihilation products of dark matter particles may be easily detectable. If photons 
are produced in dark matter annihilations, for instance, existing telescopes can detect them. Another attractive 
feature of indirect detection is that the velocity-weighted, thermally averaged annihilation cross section of the dark 
matter, (crv), which governs the expected indirect detection signal, is constrained if the dark matter is a thermal relic. 
Additionally, indirect detection has the potential to reveal information about the distribution of dark matter beyond 
our local environment. For these and other reasons, indirect detection nicely compliments the other techniques that 
may be used to identify dark matter: direct detection and collider searches @. It is likely that all three techniques 
will be necessary for a definitive identification of the dark matter. 

The present is an exciting time for indirect detection as a number of experiments are currently underway that are 
capable of detecting a signal from dark matter annihilations. Neutrino detectors such as IceCube and AMANDA H, 
air Cherenkov detectors such as H.E. S.S. and VERITAS @, cosmic ray detectors such as PAMELA [6J, and space 
based 7-ray telescopes such as the Fermi Gamma-Ray Space Telescope (FGST) [7J are all poised to make important 
contributions to the indirect detection of dark matter. 

In this paper, we focus on the possibility of using 7-ray data taken by the Large Area Telescope (LAT) on board 
the FGST to constrain the properties of dark matter. The LAT is a wide field, pair conversion 7-ray telescope that 
covers an energy range from about 20 MeV to 300 GeV [7j. Gamma-rays are particularly well suited for indirect 
detection because they are relatively easy to detect, they are produced in many models of dark matter annihilation, 
and they propagate through the universe with small optical depth (especially at low energies). 

Unfortunately, although 7-rays may in principle contain significant information about dark matter, the process of 
extracting this information is severely complicated by the presence of large and uncertain backgrounds to the dark 
matter signal. The primary challenge of indirect detection using 7-rays is therefore to extract a signal which may 
be hidden in backgrounds that are larger by orders of magnitude. Two primary features of the detected photons 
are their energy and arrival direction. A number of studies have used these two features to extract the dark matter 
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signal from the background. For example, Ref. @ used only energy information, while Ref. @ and Ref. 10 used 
both the spectral and angular distribution information. Given the photon counts, derived quantities can also be used 
to distinguish the signal from the backgrounds. The probability distribution function (PDF) has been proposed as a 
discriminant by a number of groups [Ill4l3|. Anisotropy of the distribution, especially when combined with spectral 
information, has also been proposed as a powerful way of extracting the signal [HI, Gil • 

Here we propose a new and robust approach for constraining dark matter from 7-ray data that uses only the 
angular distribution of the photons. We rely on two important aspects of the dark matter signal to help separate it 
from backgrounds. First, part of the expected signal from dark matter in our galaxy is smooth (i.e. variations from 
the mean flux in nearby pixels are uncorrelated) . Second, the dark matter signal comes from a nearly spherically 
symmetric halo, so the signal is azimuthally symmetric about the axis connecting us to the center of our Galaxy. This 
is in sharp contrast to the backgrounds, which are heavily concentrated near the disk of the Galaxy. More generally, 
astrophysical backgrounds have different morphologies and may be clumped (i.e. variations from the mean flux may 
be correlated in nearby pixels). These differences between the signal and the backgrounds allow us to remove some 
of the contribution from the backgrounds to place an interesting limit on the dark matter. The approach, which we 
call a Ring Analysis, makes very minimal assumptions about the nature of the signal and no assumptions about the 
backgrounds. Therefore, this approach is very conservative and will lead to robust limits on the properties of the dark 
matter. 

In Sjl] we describe the way we processed the Fermi LAT data to generate photon count and exposure maps. In 
£|IIII we present the Ring Analysis technique that we have developed for constraining the presence of an azimuthally 
symmetric signal on the sky. Consider an annulus centered on the axis connecting us to the Galactic Center, identified 
by the angle ip between this axis and the annulus. The Ring Analysis results in an upper limit on any contribution to 
the flux that is uniform in this annulus. Fig. U presents these upper limits on the uniform flux from the Fermi LAT 
data. Transforming these upper limits into constraints on the properties of dark matter particles requires a number 
of steps and assumptions. These, and in particular the uncertainties involved, are discussed in i jlVl Our conclusions 
are presented in JV] 



II. DATA 



The Fermi LAT is a pair conversion 7-ray detector that operates roughly in the energy range from 20 MeV to 300 
GeV. A scintillating anti-coincidence detector allows for the rejection of contaminating high energy particle events. The 
specifications of the detector are described in detail in Ref. [7| . Our analysis is based on LAT data downloaded in the 
form of weekly all-sky releases from the Fermi Science Support Center website at http://fermi.gsfc.nasa.gov/ssc/data/ 

Even with the on-board anti-coincidence detector there is a residual background of particles that are misclassificd 
as 7-rays by the LAT detectors. This poses a challenge for our analysis because these misinterpreted cosmic rays 
constitute a potentially large and uncertain background. The Fermi collaboration has made public the DataClean 
event class which implements several data cuts to minimize cosmic ray contamination as well as improved particle 
background and instrument modeling. Their data selection techniques and event modeling are described in Ref. [l6| . 
We restrict our analysis to only the events labeled as DataClean and use the corresponding P6_V3_DATACLEAN 
instrument response function to calculate exposure maps. 

Following Ref. [16j, we confine our analysis to those events coming from zenith angles < 100° in order to reduce 
contamination by 7-rays produced in cosmic ray interactions with the Earth's atmosphere. The resulting data set 
amounts to an exposure of roughly 7 x 10 10 cm 2 s across the sky covering a date range from 2008-08-04 to 2010-12-07. 

We produce photon count and exposure maps using the GaDGET package, a set of software routines designed 
for use with the LAT data by the Fermi collaboration [l2j. In generating the maps, the data were divided into 29 
energy bins between 1 GeV and 100 GeV, logarithmically spaced in energy. This restriction in energy was chosen 
because the LAT performance is well characterized in this energy range. Because the generation of these maps is 
computationally intensive, the data were processed in parallel on the Fulla 1 computing cluster at Fermilab. The maps 
generated using the GaDGET software were then converted to the HEALPix 2 isolatitude pixelization scheme with 
A^side = 64, corresponding to a pixel size of roughly (l ) 2 (comparable to the width of the point spread function 
(PSF) of Fermi at the lowest energy we consider). Fig. [1] shows the resulting all-sky 7-ray map in the energy range 
1 GeV < E < 100 GeV . 



1 http://fulla.fnal.gov/ 

2 http://healpix.jpl.nasa.gov/ 



FIG. 1: All-sky map of the 7-ray flux between 1 GeV and 100 GeV as measured by the FGST with an exposure of roughly 
7xl0 10 cm 2 s. 



III. RING ANALYSIS 
A. Overview 

The Ring Analysis technique that we develop in this paper is a method for constraining the presence of a signal 
on the sky that satisfies two requirements: 1) it is azimuthally symmetric, and 2) variations from the mean signal 
at any zenith angle are uncorrelated from one pixel to the next. As we discuss in §IV( the annihilation signal from 
smooth dark matter is expected to satisfy these two requirements, and we can quantify the extent to which the signal 
deviates from these conditions. In this section, however, we make no reference to the nature of the signal itself. 

A consequence of these two features is that the signal in each pixel of a ring of constant ip (where ip is the zenith 
angle) can be considered to be drawn independently from some underlying distribution. In the statistical literature, 
such a signal is referred to as independent and identically distributed, or i.i.d. The constraint will be on the amplitude 
of any i.i.d. component of the data. In the case of interest, the backgrounds dwarf the signal and the background flux 
in each pixel may depend on nearby pixels and may vary as a function of the angle transverse to ip, which we label <p. 

Consider a ring on the sky of constant ip defined by ipi — Aipi/2 < ip < ipi + Aipi/2, where ipi labels the central ip 
of the ith ring, and Aipi is its angular width. Since we have divided the sky into pixels, we define Fi((pj) to be the 
flux in the jth pixel (labeled by its azimuth angle <pj) of ring i. We assume that the pixels have been sorted in terms 
of <pj so that, for example, pixels j and j + 1 appear adjacent to each other on the sky; such sorting preserves any 
non-i.i.d. component of the data. The flux Fi((pj) receives contributions from an i.i.d. component due to the signal, 
F it g(4>j), and contributions from some possibly non-i.i.d. components due to the backgrounds, F it B((pj) ( see upper 
panel of Fig. [5] for an illustration). Our goal is to obtain an upper limit on the mean contribution from the signal. 

Consider the fluxes in all pixels in the ring shown in Fig. [2] This sequence is clearly not i.i.d. because of the 
peaks at (p ± 90°. That is, the flux in the pixel at cp = 90° is clearly correlated with the flux in nearby pixels. The 
physical reason for this is clear as well: these values of <p correspond to the plane of the galaxy where backgrounds are 
particularly large. The mean flux of any i.i.d. signal in this ring is clearly much smaller than the flux at the peaks. 
Choosing as a constraint the mean flux in all pixels is also not optimal as the plane of the Galaxy is distorting the 
mean. Rather, we expect the constraint on the mean i.i.d. flux to be at the level of the flux away from these peaks. 
One way to arrive at this systematically is to: 

• Start with the observed distribution (Fi((pj)) and note that it is not an i.i.d. sequence 

• Create a flux threshold (Ff) and remove all pixels with flux above the threshold 

• Test and see if the resulting (truncated) sequence (Gi((pj, F?)) is consistent with being i.i.d. 
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FIG. 2: An illustration of the terms defined in the text for the case of the data from the Fermi telescope. Top panel: The flux 
in a ring of constant ipi, Fi(cf>j), as a function of </>, the angle transverse to ip. For the purposes of illustration, we have chosen a 
ring centered at tpi = 15° of width 6°. The energy of the photons has been restricted to 1 GeV < E < 100 GeV . The increase 
in flux at (f> « ±90° is due to the galactic plane. The dashed line illustrates the maximum flux that is consistent with the the 
data being i.i.d. as determined by the BDS test. Bottom Panel: The function Gi{4>j, F^) obtained by dropping all pixels with 
flux above the threshold Fj" in the top panel. This new set of fluxes is inconsistent with an i.i.d. sequence; a slightly lower 
value of Fj would produce a Gi that is consistent with an i.i.d. sequence. Therefore, the chosen value of F T is the strongest 
constraint on the i.i.d. flux in this ring. 



• If the truncated sequence is not consistent with an i.i.d. distribution, lower the threshold and repeat 

• Once the truncated distribution is consistent with i.i.d., any i.i.d. component will generally have mean flux that 
is below the median of the truncated sequence, — median(Gi (</>.,■, Ff)) : so this median value becomes the 
upper limit on the i.i.d. flux in the ring (there are exceptions to this statement, an issue which we address 
below) 

To determine whether a given Gi(<pj, Ff) is i.i.d., we use the Brock, Dechert and Scheinkman (BDS) statistic fl8| . 
The BDS statistic tests the null hypothesis that a sequence is i.i.d. by measuring the degree of spatial correlation 
in the sequence. In essence, this is accomplished by searching for sub-sequences of length m that are significantly 
different from other m-long sub-sequences in the data; the value of m is referred to as the 'embedding dimension'. The 
null hypothesis can be rejected if the BDS statistic falls outside of some desired confidence interval; in our analysis we 
use a 3cr confidence limit. See Ref.[l9| for an introduction to the BDS statistic. To implement the BDS test we use 
a code made available in Ref. (20j . The test itself depends on two parameters: the maximum embedding dimension, 
mmaxi and a parameter e which effects how the discrepancy between different m-sequences is measured. Following 
the recommendations of Ref. [2l[, we use e = 0.5a where a is the standard deviation of Gi(4>j,F^) and m max = 5; 
we find, however, that our results are not very sensitive to the choices of these parameters. The results of the BDS 
test can be easily checked by eye; any non-i.i.d. behavior is obvious to the eye as spatial clumping. 

Ideally, the limit we derive through this method will be significantly lower than the mean flux in the ring, Fi(tfij), 
because we have effectively removed the contributions to the ring that are not i.i.d. The upper panel of Fig. [2] shows 
the initial Fi(cj)j) for a particular ring on the sky; the lower panel shows G^j^Ff), where Ff has been determined 
using the BDS test. 
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B. Monte Carlo Testing 

There are several important qualifications to the above discussion. First, it is possible to engineer pathological 
signals and backgrounds such that the limit determined by the Ring Analysis method is actually lower than the mean 
signal flux. Any realistic astrophysical sources are unlikely to have such pathological distributions, however. Second, 
the point spread function (PSF) of the Fermi telescope at 1 GeV is comparable to the size of our pixels; one might 
worry that the non-zero PSF could lead to correlations between pixels that would invalidate the i.i.d. property of 
the signal. Third, the statistical literature recommends that the BDS test be used only on data sets that are large, 
preferably with more than 500 elements [211 ] . In our analysis, however, we apply the test to data sets with as few as 
50 elements, clearly pushing the limits of the BDS test. Finally, if the signal is a smoothly varying function of ip then 
it will not have a constant value in any ring of finite angular width since different pixels at fixed ipi integrate over the 
flux with different weighting over ip. All of the above issues can be addressed through the application of Monte Carlo 
tests. By applying the Ring Analysis technique to simulated data sets, we can determine to what extent and under 
what conditions our constraints are valid. 

Since we aim only to place an upper limit on the mean flux of an i.i.d. component in the data, the relevant statistic 
for evaluating the success of the Ring Analysis technique is the ratio of the calculated upper limit on the i.i.d. flux, 
F UL , to the mean i.i.d. flux, /ii.i.d.. Ideally, i 7 ' UL / pti.i.d. will always be greater than one (so that our constraint is 
valid) but not much greater than one (so that the constraint is as tight as possible). We wish to characterize how this 
statistic varies as a function of the mean i.i.d. flux. 

We generate mock data sets by combining a simulated i.i.d. signal and a simulated background. These mock data 
sets are then analyzed using the Ring Analysis technique and its performance is evaluated. The mock i.i.d. component 
is drawn randomly from a Poisson distribution with a desired mean. The random draws ensure that the signal is in 
fact i.i.d. and the assumption of a Poisson distribution is justified for most astrophysical sources. For the background 
model, it makes sense to use the observed data itself since the data is likely background dominated. However, since 
any i.i.d. component to the background will increase F / fi- u i.d.> we subtract the determined i.i.d. flux from the 
observed background to generate our background model. This is the most conservative test of our analysis: if the true 
background contains any i.i.d. component then the true -F UL //ii.i.d. can only be larger than the -F //ii.i.d, measured 
in the Monte Carlo trials. 

One might worry that the point spread function (PSF) of the FGST could lead to correlations between pixels that 
might disturb the i.i.d. nature of the underlying signal. The width of the FGST PSF is a declining function of energy; 
at the lowest energy that we consider (1 GeV), the 68% containment angle is slightly less than 1 degree (for normal 
incidence; the PSF broadens slightly when the incidence angle is beyond about 50°). Consequently, the width of 
the PSF at the lowest energies considered is comparable to the size of our ~ 1° pixels. To account for effects of the 
PSF we applying a Gaussian smoothing kernel to the mock signal with a standard deviation of 0.5 pixels. Since the 
PSF actually gets narrower at higher energies, the application of such a smoothing kernel to all of the photon data is 
conservative; at high energies we are overestimating the size of the PSF. 

As mentioned above, we do not expect the signal in an angular ring of finite width to be exactly constant even if 
the signal is azimuthally symmetric since the signal is a smoothly varying function. To account for this possibility in 
our Monte Carlo tests we introduce a variation in the mean signal flux that amounts to a 40% decrease across the 
ring. This value is chosen because it is characteristic of the variation in the smooth dark matter signal in which we 
are ultimately interested. 

Fig. |3] presents the results of our Monte Carlo tests of the Ring Analysis method. The figure shows the quantity 
-F UL //ii.i.d. for the mock data sets as a function of /Ui.i.d./MsG; where /j,bg is the mean background flux (the signal 
and background have been simulated as described above with the effects of the PSF and the smoothly varying signal 
included). For each value of \ibg we drew 1000 realizations of the i.i.d. signal. The shaded region indicates the 
95% confidence interval for these 1000 mock data sets. Fig. [3] reveals that the Ring Analysis technique is performing 
essentially as hoped: F L //ii.i.d. is always larger than one and it stays close to one over a fairly large range of /ii.i.d.- 
For very low /ii.i.d. the signal makes essentially no contribution to the data and the flux threshold therefore stays 
constant; consequently -F UL //ii.i.d. oc l//i;.i.d. in this regime. 

C. Constraints on the i.i.d. Flux 

Fig. 2] shows the upper limits on an i.i.d. component in the Fermi data derived using the Ring Analysis technique 
as a function of the angle tp from the galactic center. As one moves away from the Galactic Center (i.e. towards 
tp = 90°), the limit becomes tighter simply because the flux is lower; beyond ip = 90° the limit becomes weaker again 
as we look toward the galactic anticenter. In generating Fig. |4]we divided the sky into 45 rings of equal width in if). 
The BDS threshold curve is shown only for those rings for which the test could be conducted on more than 50 pixels 
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FIG. 3: Results of the Monte Carlo trials for evaluating the performance of the Ring Analysis. The y-axis represents the ratio 
of the upper limit flux, F , determined using the ring analysis to the mean flux of the mock i.i.d. signal, /iu.d; the x-axis 
represents the ratio of the mean flux of the mock i.i.d. signal to the mean flux of the mock background. The shaded band 
indicates the 95% confidence interval of the Monte Carlo trials; 1000 trials were conducted for each value of /iu.d.. From the 
figure it is clear that the Ring Analysis is performing as hoped: -F UL //ii.i.d. is always larger than one (so that the constraint is 
valid) and it stays close to one over a fairly large range of ^u.d. (so that the constraint is tight). 



in accordance with the results of our Monte Carlo simulations. As can be seen from the figure, the BDS threshold is 
significantly more constraining than the mean flux in the ring. 

As a consistency check on our results, we compare the minimum threshold flux in Fig. 0] (which occurs near tj) = 90° 
as expected since this high-latitude region is least contaminated by the Galax y) t o the value of the extragalactic 7-ray 
background determined by the analysis of the Fermi Collaboration in Ref. [16[ (see the dash-dotted line in Fig.@|. 
Since any isotropic, smooth background is necessarily i.i.d. we expect our minimum threshold flux to be at least as 
large as the measured isotropic background. We do not expect the minimum threshold flux to be much larger than the 
isotropic background, however, because the isotropic background dominates at if) ~ 90°. As expected, our minimum 
threshold flux of roughly 5 x 10 _7 cm _2 s _1 sr~ 1 is slightly larger than the value of ~ 4 x 10 _7 cm~ 2 s _1 sr _1 obtained 
by Ref. [l6[ for the extragalactic 7-ray background flux over the same energy range (1 GeV < E < 100 GeV ). 

Also shown in Fig. U is the flux predicted from the smooth dark matter component for a particul ar p article mass 
and cross section (/wimp is defined in the next section) and a smooth Navarro-Frenk- White (NFW) [24] distribution 
with canonical values for the total halo mass and concentration. The shaded region (described in the next section) 
reflects the uncertainties in the dark matter distribution. Roughly, then, these values of the mass and cross section 
are ruled out by the Ring Analysis. In the next section, we will project this constraint on to the mass-cross section 
plane and propagate the uncertainties in the dark matter distribution to this plane. Here, we note that the most 
stringent limit comes from the ring 3 with ip = 15°; this model-independent limit is: 

Fu.d.O = 15°) < 4.5 x lO^cm-V^r -1 . (1) 

The fact that this analysis identifies the annulus at ip — 15° as the most constraining is consistent with the arguments 
of Ref. 



3 Note that our analysis is restricted to if> > 7°. While the constraint could in principle be improved by moving to smaller ip, the small 
number of pixels at such ip means that the BDS test loses much of its statistical power. We have confirmed with Monte Carlo tests 
that the constraints placed by the Ring Analysis tests become invalid in this regime. For the purposes of constraining the dark matter 
signal, Fig. f4] suggests that the innermost region will not tighten the constraints. 
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FIG. 4: Plot of the upper limit flux, 7^ , in the ith ring (solid line) determined using the BDS test as a function of the angle 
from the Galactic Center, ip, at which the ring lies. This threshold flux represents the maximum flux possible for any i.i.d. 
component in that ring. The dotted curve shows the mean flux in the ith ring for comparison. The dashed curve represents the 
expected signal from the smooth dark matter in our galaxy assuming /wimp = 8 x 10 -28 cm 3 s -1 GeV -2 and a canonical halo 
model; the shaded region represents the errors on the expected dark matter signal when uncertainties in the halo properties 
arc taken into account. The dash-dotted line corresponds to the value of the extragalactic flux determined by the analysis of 
Ref. We have restricted the data in the plot to 1 GeV < E < 100 GeV and have used 45 rings of equal angular width in 
ip across the sky. 



Fig. [4] also shows the mean flux in each ring. The mean flux is always larger than the i.i.d. upper limit; in the most 
constraining ring at ip — 15°, the i.i.d. limit is a factor of 2 below the mean flux. This factor of two illustrates the 
power of the i.i.d. analysis: the limit obtained from the mean in the ring is contaminated by flux near the Galactic 
plane, contamination that is removed by the Ring Analysis introduced here. 

The main advantage of the Ring Analysis is that it makes no assumptions about the source or properties of the 
backgrounds to the dark matter signal. This is a significant advantage because uncertainties in backgrounds are the 
main limitation for constraining dark matter with 7-ray observations. We have assumed only that the signal in each 
pixel is drawn from a distribution that is invariant under rotations around the Galactic Center. Of course, the cost of 
assuming little is that the limit that can be placed using this technique is comparatively weak. For instance, since any 
isotropic backgrounds present in the data will also meet these assumptions, their presence will make our determined 
limit worse. Furthermore, we note that this technique is not well suited for actually detecting dark matter, but rather 
should be viewed as a way of placing upper limits on the dark matter signal. 



IV. CONSTRAINTS ON DARK MATTER 



A. The Dark Matter Signal 



The Milky Way is believed to exist within a roughly spherical halo of dark matter [24| . Consequently, the annihi- 
lation signal from galactic dark matter is expected to have azimuthal symmetry around the line connecting us to the 
Galactic Center. The dark matter constituting the halo is in turn thought to exist in two forms: a smooth component 
and a clumped component termed subhalos |25f . Almost by definition, a smooth dark matter component will have 
an annihilation signal for which variations from the mean are uncorrelated. Therefore, the annihilation signal from 
smooth dark matter is a perfect candidate for the Ring Analysis developed in the previous section. Here we ignore 
potential signals from annihilations occurring in subhalos, as well as a possible cosmological signal from dark matter 
annihilations. Therefore, the constraint - on only the smooth Galactic component - is conservative. 



8 



In order to turn the flux limit derived by the Ring Analysis into constraints on the dark matter particle properties, 
we must develop a model for the expected annihilation flux from smooth Galactic dark matter. The observed flux of 
photons due to such annihilations along a given line of sight is 



where the first factor on the right 



FN>) = /wiMpJW (2) 

47T 



_ N,{av) 

/wimp = 2 (3) 

X 

depends on the particle physics properties of the dark matter: mass M x , velocity- weighted annihilation cross section 
(crv), and photon counts per annihilation 

N ^L,^ iE (4) 

within some desired band E m i n < E < E max . The second factor in Eq. ^ 

J(ip) ee / dldn P 2 (l,ip) (5) 



depends on the distribution of the dark matter halo, with p(l, ip) the smooth dark matter density at line of sight 
distance I and zenith angle ip. We have assumed that p is spherically symmetric about the Galactic Center so that 
the observed flux can be written as a function of ip only but will revisit this assumption later. A given profile and 
therefore a value of J translates the constraints on the i.i.d. flux in the previous section into a constraint on /wimp ■ 
A canonical J emerges from assuming an NFW profile with a scale radius of 20 kpc and local dark matter density 
of po — 0.43 GeVcm~ 3 . This leads to the J(V0 depicted as the dashed line in Fig. 2] for the given value of /wimp • 



B. Uncertainties in the Dark Matter Profile 



The dark matter profile is not constrained very tightly by observations, and this uncertainty propagates to the 
constraints on the particle physics properties. We address this in two ways here, the first applies to all constraints on 
the smooth halo and the second is specific to our assumption that the dark matter distribution is spherical. 

To allow for freedom in the dark matter profile, we piggy-back on the analysis of Ref. [26j. They used nine sets of 
observational data to constrain the dark matter profile in our Galaxy. They assumed a density of the form 



.2, 



o(r) = - — - f6) 

M ' AnalG (r/a h )(l+r/a h r-y 



where r is the distance from the Galactic Center, ah is the scale radius of the dark matter halo, ah is a scale velocity, 
and 7 is a parameter which effects the behavior of the central density cusp, with 7=1 corresponding to the usual 
NFW form. A given set of these free parameters translates into a set of values for J(ip) assuming the distance between 
the Earth and the Galactic Center to be 8.5 kpc. 

The allowed values of the three profile parameters (and therefore J(ip)) are contained in Monte Carlo Markov 
Chains run in Ref. [26]. They have kindly provided us with those chains, and Fig. [5] projects these on to J(ip) for 
three different values of ip. It is clear from this figure that the width of the distributions increases rapidly with 
decreasing ip. 

The widths of the distributions show in Fig. [5] can be propagated directly into error bars on the expected signal 
for dark matter in Fig. @] The shaded region in Fig. @] represents the 95% confidence interval for the allowed dark 
matter signal for a particular choice of /wimp ■ The value we have chosen for /wimp is illustrative in the sense that 
it is roughly the lowest value that is excluded by the Ring Analysis limit. 

Our analysis so far has assumed that the signal from the smooth component of the dark matter is spherically 
distributed around the Galactic Center. In fact, the true shape of the Milky Way's dark matter halo may be triaxial 
27]. If this is the case, then the dark matter signal is not uniform in rings of constant ip on the sky; instead, the signal 
in such a ring will be an oscillating function of the azimuth angle 4> (the thick curves in Fig. [5]). As a result, the total 
signal in a given ring will no longer be i.i.d. (i.e., it is not identical). However, there will still be some component of 
the signal that is i.i.d.; the flux of this component is given by the minimum flux of the dark matter signal in the ring 
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FIG. 5: Histogram of values of J(f/>) = J dlp 2 {l,ip) corresponding to models in the MCMC chains provided by [26l ]. We show 
the histograms for ip = 10° and ip = 15° since these angles provide the strongest constraint on dark matter annihilations. The 
broad histogram for ■0 = 5° illustrates the well-known statement that the uncertainty in the flux increases towards the Galactic 
Center. 



(the thin horizontal lines in Fig. [5]). Since the triaxiality of the halo effectively decreases the magnitude of the i.i.d. 
signal, the stated lower bound is too aggressive. It is therefore important to quantify how the triaxiality of the halo 
effects the level of the i.i.d. signal. 

In order to estimate the magnitude of this effect, we calculate 

J(ip,cf>) = J dldflp 2 (l,^ : (j)) (7) 

along different lines of sight in model triaxial galaxies. J(ip, (j>) is the analogue of JW) for a dark matter halo that is 
not assumed to be spherically symmetric. Observations in the Milky Way by Ref. |27J and simulations of disk galaxies 
by Ref. [28| suggest that one of the axes of the inner halo (which is also the region that dominates our exclusion limit) 
is aligned with the rotation axis of the galactic disk. Following the results of Ref. 27], we assume that the minor 
and major axes of the halo lie in the galactic plane; we find, however, that this choice of alignment does not have a 
significant impact on our results. We assume a density profile of the form Eq. [5] with the replacement 

t 2 = (x/cf + (y/af + {z/bf (8) 

where the x, y and z axis are aligned with the minor, major and intermediate axes of the halo respectively. Strictly 
speaking, Refs. [29], [3(| have found that the axis ratios a/b and b/c do not remain constant throughout the halo; 
rather, they decrease slightly (i.e. the ellipsoid becomes more elongated) towards the center of the halo. However, 
since our exclusion limit is dominated by a small range of distances from the Galactic Center we feel justified in 
approximating the halo by ellipsoids with constant axis ratios and alignments throughout the galaxy as per Eq. [5] 

Fig.[S]shows J(ip — 10°, 4>) along different lines of sight in model galaxies with several values of the axis ratios (thick 
lines). In generating this plot we have used ah = 270km/s, at = 5.9 kpc and 7 = 0.028; in a spherical halo these 
values yield a J(ip — 10°) which is at the fifth percentile of all those calculated from the Markov chains of Ref. [26| . 
Since the uncertainty in our exclusion limit is dominated by uncertainty in the halo model, this parameter choice is 
therefore appropriate for our 95% exclusion limit on the dark matter properties. We plot only those lines of sight with 
■ijj = 10° because this is the angular range relevant to our exclusion limit. All of the curves have been normalized so 
that they have the same mean as the corresponding curve for a spherical halo; since the total amount of dark matter 
in the inner halo is strongly constrained (by measurements of galactic rotation curves, for instance), we expect this 
choice of normalization to be reasonably accurate. 
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To proceed further, we need to know the axis ratios of the isodensity surfaces throughout the Milky Way halo. 
Ref. [27| determined values for the axis ratios of some of the isovelocity surfaces in the Milky Way. While these values 
could in principle be converted to axis ratios for the isodensity surfaces given a density model, the large uncertainties 
associated with these measurements lead us to take a hopefully more robust approach. Using numerical simulations, 
Refs. [29j, |30j have determined probability distributions for the axis ratios of dark matter halos. We run Monte Carlo 
simulations based on these results in order to quantify the extent of the variation in the dark matter signal due to 
triaxiality. 

We calculate J(ip, 4>) for halos with axis ratios drawn from the probability distributions of Eqs. 17 and 18 in 
Ref. [Hj]. We adjust these distributions slightly to account for the fact that they are calculated at a distance from the 
Galactic Center, r 2 5oo, given by pfosocO/Pcrit = 2500 (note that r25oo is defined in terms of the local halo density and 
not the mean interior density) while we are interested in the distributions much closer to the Galactic Center. Using 
Eqs. 6 and 7 of Ref. [2!| we scale the axis ratios so that they correspond to a distance of ~ 1.5 kpc from the Galactic 
Center, roughly the innermost distance that effects our constraint. This corresponds to decreasing a/c and b/c by 
roughly 30% and 10% respectively. The resultant distributions are roughly consistent with the results of Ref. [30j j . 
Finally, we assume that the line connecting the sun to the Galactic Center lies in the plane of the minor and major 
axes of the halo but with a random orientation in that plane. 
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FIG. 6: The value of J(l/i = 10°,</>) along different lines of sight for three model triaxial galaxies (thick lines). The minor, 
intermediate and major axes are labelled by a, b and c respectively. Also plotted are the corresponding levels of the i.i.d. 
component of the signal (thin lines) . Following the results of Ref. [23] we have assumed that the minor and major axes of the 
halo lie in the galactic plane and that the line connecting the sun to the Galactic Center is offset from the minor axis of the 
halo by 15°. 



Fig. [5] shows the flux as a function of azimuthal angle ^ in a given ring for three different sets of (a,b,c). The 
key take-away from these plots is the difference between the mean flux (i.e. the mean value of J(tjj = 10°, <p) for 
the thick curves in Fig. [5]) and the minimum flux (i.e. the thin curves in Fig. [S]) in the ring. This difference is the 
amount by which we have been implicitly overestimating the i.i.d. signal. Dividing the difference by the mean flux 
then provides an estimate of the relative over-estimation, the amount by which we should loosen the i.i.d. constraint. 
Fig. [7] shows the distribution of this fractional difference (also at ijj = 10°) for 10000 Monte Carlo realizations of the 
halo axis ratios drawn from the probability distributions of Ref. [29j . It is clear from this figure that error introduced 
by assuming a spherical halo is typically no more than 25%. Since the uncertainty introduced into our exclusion limit 
by the uncertainties in ah, a,h and 7 is much greater than this, and since the degree of triaxiality in our own halo is 
not very well constrained, we choose to ignore this source of uncertainty in our exclusion plots. 

We conclude this section by mentioning that the presence of baryons may have a non-negligible impact on the shape 
of the dark matter halo, particularly in the innermost region [3lj . The effects of baryons have not been included in 
the simulations of Ref. [29] nor Ref. 30] . While there may be a disturbance due to baryons near disk, the effect of 
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FIG. 7: Distribution of the percentage variation (mean to minimum) of the signal from smooth dark matter at ip = 10° in 
10000 realizations of our Monte Carlo modeling of triaxial halos. Axis ratios have been drawn from the probability distributions 
of Ref. 29] (scaled according to the prescription described in the text). The line connecting the sun to the Galactic Center is 
assumed to lie in the plane of the minor and major axes of the halo, but with a random orientation in that plane. 

such a disturbance on our exclusion limit is likely negligible since the limit is dominated by the region away from the 
disk where the backgrounds are lowest. 

C. Constraints on Dark Matter 

The previous discussion has constrained /wimp over the energy range lGeVto 100 GeV. We find that our final 
constraint over this energy range is 

/wimp < 5.8 x l(T 28 cm 3 s" 1 GeV" 2 (9) 

In order to turn the limits on /wimp into constraints on the dark matter particle properties we must assume an 
annihilation channel for the dark matter so that JV 7 (Eq. 2]) can be calculated. The constraint on /wimp and the value 
of _/V 7 can then be transformed into a constraint in the M x -(av) parameter space using Eq. [3] 

The choice of the energy range that we impose on our analysis can have a significant impact on the constraints 
that we place in the M x -(av) plane. For E <C M x the annihilation photon spectrum is relatively flat compared to 
that of the backgrounds, so increasing E m \ n in this regime effectively increases the size of the signal relative to the 
total flux, thereby improving the constraint on the dark matter properties. As E approaches M x , the spectrum of 
annihilation photons falls off very quickly so increasing i? m i n in this regime causes the constraint to become very 
weak. We expect, then, that the optimal constraint on the dark matter comes from choosing E m i n to be some fixed 
fraction of M x . Rather than enforcing the optimal E m { n for each M x exactly (which would require dividing the data 
into many more energy bins than the 29 that we employ), we allow E m { n to vary freely for each M x that we consider. 
We then choose the value of E m i n that maximizes our constraint on the dark matter. Random fluctuations in the 
strength of the constraint induced by varying i? m i n are small, so choosing E m i n to maximize the constraint does not 
decrease the statistical significance of our result. Varying i? m i n is our one minimal use of the energy information for 
the photons; by making our analysis essentially independent of detailed spectral information we make our constraints 
more robust. 

Fig. [3] shows how the constraints that we place on /wimp in different energy bins with the Ring Analysis are 
translated into constraints in the M x -(av) plane for different assumed annihilation channels (each of which correspond 
to a different iV 7 ). We consider two possible annihilation channels: XX ~^ bb and \X ~~ ^ t t, where x is a neutralino. 
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Nj is calculated for these different channels over a mass range of 10 GeV to 10 TeV using the DarkSUSY 4 package [32| . 
Since our technique places a limit on /wimp = N 7 {av) /M^ we expect the limit on (av) to go roughly as M^. The 
flattening of the limit at low mass is due to a decrease in Y 7 as more and more annihilation photons fall outside of 
the energy limits of our analysis. 

The dashed curves in Fig.|8]show what the limit would be if the dark matter profile, as quantified by J(i/>), produced 
the mean annihilation signal, and there was no uncertainty in the profile. The uncertainty in dark matter profile then 
loosens the constraint on the annihilation cross section by more than a factor of 2. 




10" 10 s 10" 10 7 10" 10 s 10 6 10 7 

Mass (MeV) Mass (MeV) 

FIG. 8: Exclusion plots in the (<jv)-M x plane generated by the Ring Analysis for different annihilation channels of the dark 
matter. The solid region represent the 95% exclusion range with respect to the uncertainties in the properties of the smooth 
halo. The dashed line represents the boundary of the excluded region when the mean value of J{ip) from the Monte Carlo 
chains of Ref. [2(| is used to calculate the expected dark matter signal. The left panel shows the exclusion limit for neutralinos, 
X, that annihilated to bb; the right panel corresponds to neutralinos annihilating to rf . 

The limits presented in Fig. [5] are comparable to the corresponding plots from Ref. @, which also considers the 
annihilation signal from the smooth galactic dark matter component, but which makes a stronger set of assumptions 
about the signal and backgrounds. Our limits are complementary in the sense that they are almost entirely independent 
of any assumptions about the diffuse 7-ray background. We also find that our constraints are comparable to those 
obtained from an analysis of extragalactic dark matter annihilations by Ref. Q. As evidenced in Fig. 5 of that work, 
the constraints derived from extragalactic annihilations are subject to very large uncertainties in the magnitude of the 
dark matter signal. While uncertainties in the signal from smooth galactic dark matter are also important, they are 
not nearly as large. Finally, our limits are also comparable to (and competitive with) those derived from an analysis 
of the annihilation signal from dwarf spheroidal galaxies by Ref. [331 ] . 

V. SUMMARY 

We have developed a technique for constraining the presence of a smooth, azimuthally symmetric signal on the sky. 
We showed with Monte Carlo simulations that the technique is robust (i.e. the limits derived from the technique are 
never below the mean signal flux). By applying the technique to data from the Fermi Gamma- Ray Space Telescope we 
derived a constraint (Eq. (fT])) on the presence of any smooth 7-ray signal that is symmetric with respect to rotations 
about the axis connecting us to the Galactic Center. When combined with a model for the signal due to annihilations 
of the smooth dark matter component in our galaxy, these limits allow us to place constraints on the dark matter 
particle mass and cross section. While our limits are slightly weaker than other recent results, they have the advantage 
of making essentially no assumptions about the backgrounds to the dark matter signal. This is a significant advantage 
because the uncertainties in models of 7-ray backgrounds are large and often unknown. 

There are several ways these limits can be improved. Tighter constraints on the dark matter profile would reduce 
the uncertainty in J(ip), which currently degrades the ultimate limit by more than a factor of 2. Including other 
sources of signal, in particular the contribution from extra-galactic halos, would improve the signal, but in most 
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models, this contribution is smaller than that from the Galactic halo. A deeper understanding of the backgrounds 
could also be used in conjunction with this method to push the limits down further. 
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